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The basic quantity for the description of the statistical properties of physical systems is the density 
of states or equivalently the microcanonical entropy. Macroscopic quantities of a system in equi- 
librium can be computed directly from the entropy. Response functions such as the susceptibility 
are for example related to the curvature of the entropy surface. Interestingly, physical quantities 
in the microcanonical ensemble show characteristic properties of phase transitions already in finite 
systems. In this paper we investigate these characteristics for finite Ising systems. The singularities 
in microcanonical quantities which announce a continuous phase transition in the infinite system are 
characterised by classical critical exponents. Estimates of the non-classical exponents which emerge 
only in the thermodynamic limit can nevertheless be obtained by analyzing effective exponents or 
by applying a microcanonical finite-size scaling theory. This is explicitly demonstrated for two- and 
three-dimensional Ising systems. 



1 Introduction 

In the study of the static properties of finite systems the main quantity of 
interest is the density of states £1 as it contains the complete physical informa- 
tion of the investigated system pQ. To be specific, consider a magnetic system 
where f2 = Q(E, M) is a function of the energy E and of the magnetization M 
(which can be a scalar or a vector quantity, depending on the system). With 
the density of states at hand one may proceed in the canonical way by com- 
puting the free energy from which other quantities such as the susceptibility 
of the specific heat can be derived. In the last years, however, an increasing 
number of groups have chosen to directly investigate the density of states or, 
equivalently, the microcanonical entropy S 1 = lnfi (setting ks = 1)- The origin 
of this increasing interest in the microcanonical approach can be found in the 
observation that ensembles are not always equivalent 0HI3HIHliniEllEllHl • Tne 
inequivalenve of ensemble is for example encountered in situations where the 
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macroscopic limit (for example due to the presence of long-range interactions 
such as the gravitational force) does not exist. But also in finite systems with 
short-range interactions the canonical and the microcanonical ensembles may 
give different results although the ensembles have to become equivalent in the 
infinite volume limit. This is especially the case for finite systems which in the 
macroscopic limit undergo a phase transition. It is the purpose of this paper 
to discuss the finite-size signatures of phase transitions in the microcanonical 

The paper is organized in the following way. In Section 2 we recall that in the 
microcanonical approach all quantities of interest can exclusively be expressed 
as derivatives of the entropy [I]. In Section 3 we briefly discuss the case of 
discontinuous phase transitions which lead to a characteristic behaviour of 
the specific heat in finite microcanonical systems. Section 4 is the main part 
of this paper. It is devoted to the signatures of continuous phase transitions 
encountered in the microcanonical analysis of small systems. We also discuss 
different ways for extracting the values of critical exponents directly from the 
derivatives of the microcanonical entropy. Conclusions are given in Section 5. 



2 Microcanonically defined quantities 

We use in the following the language of magnetism and consider a magnetic 
system with N degrees of freedom isolated from any environment. The state 
of the system is then characterized by its energy E and its magnetization M. 

In the microcanonical analysis one considers the microcanonical entropy 
density 



as function of the energy density e = E/N and the magnetization density 
m = M/N . The entropy surface s(e, m) still exhibits a weak size dependence 
which is not explicitly denoted here. Sometimes one is only interested in the 
energy dependence and therefore studies the reduced specific entropy 
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Microcanonical quantities can be expressed by functions of derivatives of the 
specific microcanonical entropy. For example, the microcanonical temperature 
reads 

T-(^>y\ ,4, 
which leads to the following expression for the specific heat: 

' ds R (e) s 



de 2 

It is worth noting that the microcanonical temperature and specific heat can 
also be defined through the entropy density (|T|). These differently defined quan- 
tities get identical in the infinite volume limit. We refer the reader to |15j for 
a thorough discussion of this point. 

Further quantities discussed in the following include the spontaneous mag- 
netization and the susceptibility. The spontaneous magnetization m sp ^(e) is 
defined to be the value of m where the entropy Q at a fixed value of e has 
its maximum with respect to m |1UU121ITd1I17| : 

m sp N (e) s(e, m sp , N (e)) = maxs(e,m). (6) 

m 

To better understand this definition consider the density of states 
Q(eN,mN) = exp(A r s(e, m)) for a given energy e. The density of states ex- 
hibits a sharp maximum as the overwhelming majority of accessible states 
belongs to the value m 3P) N of m where s has its maximum. This is the rea- 
son why we identify m sP: N with the spontaneous magnetization of the finite 
system. This definition of the microcanonical spontaneous magnetization en- 
sures that we recover the canonically defined spontaneous magnetization in 
the macroscopic limit. Note also that one can define a magnetic field for the 
microcanonical ensemble which appears as the conjugate variable to m. At the 
spontaneous magnetization m sp> N the associated magnetic field becomes zero, 
see ^U] for further details. 

Of further interest is the suceptibility [TUj 

ds d 2 s 

x (e, m) = ? (7) 

de 2 dm 2 \ dechn ) 



where the denominator on the right hand side is in fact the local curvature 
of the entropy surface. As it stands, Eq. Q is only valid in cases where the 
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order parameter is a scalar. For an order parameter with vector character (as 
encountered for example in classical spin models like the XY or the Heisen- 
berg model) the susceptibility is given by a matrix [181116] . The entries of the 
susceptibility matrix are then dependent on the local curvature with respect 
to the different order parameter components. 

Closing this Section, we reemphasize that the microcanically defined quan- 
tities all fulfill the important request that they become identical to the better 
known canonical quantities in the infinite system after the trivial transfor- 
mation from the natural variables e and m of the micro canonical description 
to the natural variables temperature T and magnetic field h of the canonical 
approach. 



3 Discontinuous phase transitions 

The interest in the microcanonical analysis really started with the observation 
that discontinuous phase transitions display typical and very interesting finite- 
size- signatures in the microcanonical treatment [31411] . Indeed, a discontinuous 
phase transition reveals itself in finite systems by a convex intruder in the mi- 
crocanonical entropy which originates from states of coexistence of different 
phases. In the canonical treatment these states are unstable and cannot be 
accessed by equilibrium methods. Interestingly, this convex intruder leads to a 
characteristic backbending of the caloric curve and thus to a negative specific 
heat. The appearance of a negative specific heat may look like an oddity of the 
microcanonical approach. However, there have been recent claims that nega- 
tive specific heats had been measured (even so in a somewhat indirect way) in 
experiments on nuclear fragmentation ^H] or on the solid-liquid transition of 
atomic clusters |2u] . 



4 Continuous phase transitions 

The main intention of this paper is to discuss the typical finite-size signatures 
of continuous phase transitions encountered in the microcanonical approach. 
We thereby focus on classical spin systems and shall discuss in the following the 
two- and three-dimensional Ising models on hypercubic lattices with periodic 
boundary conditions. The Ising model is defined by the Hamiltonian 




(8) 
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where the lattice site i is characterized by the classical variable Si = ±1. 
The sum in Equation (|SJ) is over nearest neighbour pairs. The two-dimensional 
three-state Potts jl7lll4lfTH] and the three-dimensional XY |TB] models have 
also been studied recently in the present context, but we will not discuss these 
models here. All the mentioned models have in common that the corresponding 
infinite system displays a continuous phase transition separating a paramag- 
netic phase at high temperatures (energies) from a ferromagnetically ordered 
phase at low temperatures (energies). 

4.1 Signatures in small systems 

In the canonical ensemble the typical features of spontaneous symmetry 
breaking like the abrupt onset of the order parameter or a diverging sus- 
ceptibility are only encountered in the infinite system. In any finite system 
one does not observe diverging quantities but only rounded maxima. This 
is, however, completely different in the microcanonical ensemble where the 
features of spontaneous symmetry breaking turn up already for finite sys- 
tems [IHl[E2[I3[E3[E3[Rn[2I] . 

It is instructive to first inspect directly the entropy surface s(e, to). Figure ^ 
displays this entropy surface for a three-dimensional Ising model with 6 3 spins. 
The fully magnetized states at e = —3 and to = ±1 form the two degenerate 
gound states, whereas the fully disordered macrostate at e = m = has 
the highest degeneracy of all. This entropy surface has been computed by a 
very efficient numerical method based on the concept of transition variables 
|12j . More insights are gained by looking at cuts through the entropy surface 
along the magnetization at a fixed energy e, see Figure 01 One observes that 
the entropy exhibits only one maximum at to = for energies e above an 
energy e C) N, whereas for energies e < e Cj 7v two maxima are seen at nonzero 
magnetizations, separated by a minimum at to = 0. At the energy e Cj 7v the 
curvature in magnetization direction vanishes at to = 0. 

Recalling that the spontaneous magnetization is identified with the mag- 
netization where for a fixed energy the density of states is the largest, one 
can directly extract the value of the spontaneous magnetization at a given 
energy from cuts like those shown in Figure [2j The resulting behaviour of 
the spontaneous magnetization is discussed in Figure 01 for three-dimensional 
Ising models containing 6 3 and 8 3 spins. Two different regimes are observed: 
for e > e C) N the spontaneous magnetization is zero, whereas for e < e C) M it is 
given by to = ±m sPj 7v, in accordance with the appearance of two maxima of 
the entropy surface at lower energies. At the well-defined size-dependent en- 
ergy g c ,n the order parameter sets in abruptly already in finite systems. The 
energy e C) N may therefore be called transition energy. 

Another remarkable finite-size signature of continuous phase transitions is 




Figure 1. Entropy surface s(e,m) of the three-dimensional Ising model with N = 6 3 spins as 
function of the energy e and of the magnetization m. 

shown in Figure 21 A divergence of the magnetic susceptibility is observed for 
finite systems at the same energy e Cj jv at which the order parameter bifurcates. 
This divergence is in fact readily understood by recalling, see Equation (|7|. 
that the susceptibility is proportional to the inverse of the curvature of the 
entropy surface. It is the existence of a point of vanishing curvature that is 
responsible for the observed divergence of the susceptibility, see Figure EJ 

The microcanonical specific heat, on the other hand, does not diverge in 
finite systems. It nevertheless also displays an intriguing finite-size behaviour, 
as revealed recently for models with a diverging specific heat in the infinite 
system |15j . Indeed, in the microcanonical ensemble the maximum of the spe- 
cific heat of the corresponding finite systems varies non-monotonically when 
increasing the size of the system. Starting with small systems, it first decreases 
for increasing system sizes, before it increases for larger systems and finally 
diverges in the macroscopic limit. This behavious is in strong contrast to the 
treatment in the canonical ensemble where a monotonic increase of the maxi- 
mum of the specific heat with the size of the system is always observed p2U25j . 
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Figure 2. Cut through the entropy surface at various fixed values of e. For e > e c jv the entropy 
exhibits a maximum at m = 0. For e < e c j\r, two maxima at nonzero magnetizations are visible. At 
B = e c ,N (dashed line) the extremum at m = changes from a maximum at higher energies to a 
minimum at lower energies so that the curvature changes its sign. 

In a phenomenological theory has been proposed that acccounts for the 
observed non-monotonic variation of the maximum of the microcanonical spe- 
cific heat. 



4.2 Determination of critical exponents directly from the density of 
states 

At this stage one may wonder whether the typical features of spontaneous 
symmetry breaking observed in finite microcanonical systems are described 
by the true critical exponents of the infinite system. A rapid inspection of 
the numerical data reveals that this is not the case. Indeed, the onset of the 
order parameter or the divergence of the suceptibility in finite microcanonical 
systems are in all cases governed by the classical mean-field exponents |10[ 
I12U17| . This is readily understood by noting that the entropy surface can be 
expanded in a Taylor series in the vicinity of e C) N (for a discrete spin model 
a suitable continuous function that represents the data has to be chosen). 
As a consequence one ends up with a Landau-like theory |17j which for the 
microcanonical entropy of a finite system is in principle exact, thus yielding 
classical values for the critical exponents. In recent years we have developed 
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Figure 3. Spontaneous magnetization m sp jv as function of the energy e measured in 
three-dimensional Ising models with 6 3 and 8 3 spins. At the size-dependent energy e c jy a sharp 
onset of the order parameter is observed. 



different approaches that enable us nevertheless to obtain the true critical 
exponents by exclusively analyzing the density of states of finite systems \V2i 

EE]. 

Before discussing these approaches for the spontaneous magnetization in 
some detail we have first to pause briefly in order to recall that in the entropy 
formalism considered here critical exponents often differ from their thermal 
counterparts. This is the case when the specific heat diverges algebraically 
at the critical point of the infinite system, implying that the critical expo- 
nent a, governing the power-law behaviour of the canonical specific heat in 
the vicinity of the critical point, is positive. It is then easily shown that the 
critical exponents x e appearing in the microcanonical analysis are related to 
the thermal critical exponents x by x e = x/(l — a) 10 . However, when the 
specific heat does not diverge but presents a cusp-like singularity as a < 0, 
we have x = x £ |18U16| . In the case of a logarithmically diverging specific heat 
the thermal and microcanonical exponents are also identical. 

Let us now come back to the determination of the critical exponents. As 
already mentioned we obtain classical values for the critical exponents when 
expanding the finite-size quantities around the energy e c ,N- The spontaneous 
magnetization, for example, displays a root singularity for energies e < e Ct N 
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Figure 4. Magnetic susceptibility x( e > m sp,jv) as function of the energy e measured in 
three-dimensional Ising models with 6 3 and 8 3 spins. The microcanonically defined reponse 
function diverges at e Ci jv already in finite systems. 

close to e Cj Ar: 

m sP)N (e) = A N (e Cj 7v - e) 1/2 (9) 

where the amplitude An depends on the system size. In the infinite system one 
expands around the true critical point e c , leading to the asymptotic behaviour 

m sp (e) = AeP' (10) 

for small e with e = £c ~ e , where e„ denotes the ground state energy per 

degree of freedom. The critical exponent (3 £ is given by (3 £ = f3/(l — a) with j3 
being the usual thermal critical exponent characterizing the singularity of the 
spontaneous magnetization. Here we suppose that a > 0, as it is the case for 
the Ising model. 

Away from the critical point the spontaneous magnetization of the infinite 
system is not any more given by the simple expression (|lUj) , as correction terms 
get important. One way to deal with these correction terms is to analyse the 
energy dependent effective exponents 



(11) 
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that yields the true critical exponent (3 £ in the limit e — ► |12| . 

In order to use Equation we need data free of finite-size effects. On 
the other hand one observes in Figure |3] that there exists a energy interval for 
which the spontaneous magnetizations computed in systems of different sizes 
are identical. This suggest the following procedure for computing j3 £ . Cal- 
culate the spontaneous magnetization for a series of increasing system sizes 
{Nj,j = 1, ■■■,n}. In the energy range where and N i+ i yield the same 
exponent (3 £e ffi select this exponent for the infinite system, until finite-size 
effects show up. Then plot the common exponent for system sizes iVj+i and 
iVj_|_2 until these begin to disagree and so on. This procedure is rather cumber- 
some but ensures that data essentially free of finite-size effects are obtained. 

The result of this approach ^5] for the two- and the three-dimensional Ising 
models is shown in Figure 21 For the three-dimensional case we observe that 
Pe,eff rapidly approaches a constant value that perfectly agrees with the ex- 
pected asymptotic value (3 £ = 0.367. In two dimensions the agreement is not 
as perfect, even so a reasonable estimate of /3 e is obtained by linearly extrap- 
olating the existing data. The situation is indeed more complicated in the 
two-dimensional Ising model, as one is dealing here with the rather patho- 
logical case of a logarithmically diverging specific heat with a = 0. In fact, 
we expect for any model with a/0a similar good agreement between the 
expected value of j3 e and the value obtained from linearly extrapolating @ e ff,e 
as for the three-dimenisonal Ising model. This expectation has recently been 
confirmed for the three-dimensional XY model |16j . 

The determination of critical exponents via effective exponents has the major 
drawback that data free of finite-size effects have to be obtained, which means 
that systems of increasing sizes have to be simulated on approach to the critical 
point. For Figure our largest systems contained 700 x 700 spins in two and 
80 x 80 x 80 spins in three dimensions. Simulating these systems is very time 
consuming even with the best algorithms available. The second approach for 
determining the values of critical exponents discussed in the following does not 
need large systems, but instead it explicitly takes advantage of the finite-size 
effects. 

In the canonical ensemble finite-size scaling theory is a valuable and often 
used tool for extracting critical exponents from finite-size data. Its starting 
point is the observation that in the asymptotic limit L — ► oo and T — ► T c 
(T c being the critical temperature) the behaviour of finite-size quantities is 
governed by scaling functions that are basically determined by the ratio L/£(T) 
with £(T) being the correlation length and L = N l / d the linear extension 
of the d-dimensional system |23|. I n EH we have developed a microcanonical 
finite-size scaling theory that takes advantage of the existence of a well-defined 
transition point e C; 7v in finite microcanonical systems. 

The starting point is the scaling behaviour of the entropy of finite systems 
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Figure 5. Effective exponent fie eff versus reduced energy e for the two- and the 
three-dimensional Ising models. The dashed lines indicate the expected values of f} e , see text. Error 
bars are only shown when they are larger than the sizes of the symbols. 



considered as a function of the energy, the magnetization and the inverse 
system size. This then leads to the scaling form 

L^m sPjN (e C)JV - e) « W (c{e c , N - e)L 1 /^) (12) 

of the spontaneous magnetization Q. Here W is a scaling function charac- 
terizing the given universality class and C is a nonuniversal constant which 
is different for the various model systems belonging to the same universality 
class. The exponent v e is given by u E = u/(l — a) where v is the usual ther- 
mal critical exponent governing the divergence of the correlation length on 
approach to the critical point. As the spontaneous magnetization presents a 
square root singularity in any finite microcanonical system, the scaling func- 
tion W varies for small scaling variables x = C(e Cj 7v — e)L 1//l/e as W(x) ~ yfx. 
In FigurelHlwe test the scaling form (J12[) for the usual Ising model (JHJ) with only 
nearest neighbour interactions (denoted by I) and for a generalized Ising model 
with equivalent nearest and next-nearest neighbour interactions (denoted by 
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Figure 6. Microcanonical finite-size scaling plot for the Ising model with only nearest neighbour 

interactions (I) and for the Ising model with equivalent nearest and next-nearest neighbour 
interactions (Gl). Critical exponents are determined by the best data collapse. By adjusting the 
nonuniversal constant C the data of both models fall on a unique curve, thus demonstrating the 
universality of the scaling function W given in Eq. 1121 . Error bars are smaller than the sizes of the 

symbols. 

Gl). The Hamiltonian of this latter model is given by 



7~Lgi = — ^2 SiSj — SiSk 

(hj) (i,k) 



(13) 



where the second sum is over bonds connecting next-nearest neighbours. The 
model (fT3)) belongs to the same universality class as the standard Ising model 
(JHJ) so that both models should yield the same values for universal quantities. 

Figure El illustrates two different aspects of the scaling behaviour of the mi- 
crocanonically defined spontaneous magnetization. On the one hand it shows 
that both for the Ising model and for the Gl model a data collapse can be 
achieved by plotting L^^ Vf - function of (e Cj 7v — e)L x l Vc . The values of 

the involved critical exponents are thereby obtained by the best data collapse. 
This yields the values (3 £ /u £ = 0.54 ±0.03 (0.51 ±0.03) and l/u e = 1.43 ±0.01 
(1.43±0.01) for the Ising (Gl) model, in excellent agreement with the expected 
values /3 e /v e = 0.52 and l/u £ = 1.43. On the other hand Figure El also proves 
that the scaling function W appearing in Eq. (|12|) is the same for both models, 
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in accordance with its expected universality. Indeed, the data of both models 
fall on a common master curve when adjusting the nonuniversal constant C, 
so that Cqi = 0.219C/. As shown in ^3] scaling functions obtained for models 
belonging to different universality classes are indeed different. 

The microcanonical finite-size scaling theory, which has also been applied 
successfully to the Potts ^3] and to the XY ^5] models, is very promising. One 
of the remarkable point is that no a priori knowledge of the infinite system is 
needed. Especially, one does not need to know the exact location of the critical 
point of the infinite system. 



5 Conclusion 

In this paper we have discussed the signatures of phase transitions in finite 
microcanonical systems. Discontinuous phase transitions are announced in 
the microcanonical ensemble by a backbending of the caloric curve and the 
appearance of negative specific heat. Continuous phase transitions also lead 
to intriguing finite-size signatures. Indeed, typical features of spontaneous 
symmetry breaking like the onset of the order parameter or the divergence of 
the susceptibility are already encountered in finite microcanonical systems. 
These singularities are governed by the classical mean-field exponents. We 
have discussed in this paper two different approaches that nevertheless enable 
us to extract the true critical exponents from the density of states of finite 
systems. One approach involves effective exponents, whereas the other is 
based on a microcanonical finite-size scaling theory. The latter approach 
explicitly takes advantage of the existence of a well-defined transition point 
in any finite microcanonical system. 

Acknowledgements: 

It is a pleasure to thank Alfred Huller for many years of fruitful and 
enjoyable collaboration. 

References 

[1] D. H. E. Gross, Microcanonical Thermodynamics: Phase Transitions in 'Small' Systems (Lecture 

Notes in Physic 66) (World Scientific, Singapore 2001). 
[2] W. Thirring, Z. Phys. 235, 339 (1970). 
[3] M. Bixon and J. Jortner, J. Chem. Phys. 91, 1631 (1989). 
[4] A. Huller, Z. Phys. B 93, 401 (1994). 

[5] R. S. Ellis, K. Haven and B. Turkington, J. Stat. Phys. 101, 999 (2000). 
[6] T. Dauxois , P. Holdsworth and S. Ruflb, Eur. Phys. J. B 16, 659 (2000). 
[7] J. Barre, D. Mukamel and S. Ruffo, Phys. Rev. Lett. 87, 030601 (2001). 
[8] I. Ispolatov and E. G. D. Cohen, Physica A 295, 475 (2001). 
[9] F. Gulminelli and Ph. Chomaz, Phys. Rev. E 66, 046108 (2002). 



14 



[10] M. Kastner, M. Promberger and A. Huller, J. Stat. Phys. 99, 1251 (2000). 
[11] D. H. E. Gross and E. V. Votyakov, Eur. Phys. J. B 15, 115 (2000). 
[12] A. Huller and M. Pleimling, Int. J. Mod. Phys. C 13, 947 (2002). 
[13] M. Pleimling, H. Behringer and A. Huller, Phys. Lett. A 328, 432 (2004). 
[14] J. Hove, Phys. Rev. E 70, 056707 (2004). 

[15] H. Behringer, M. Pleimling and A. Huller, J. Phys. A: Math. Gen. 38, 973 (2005). 
[16] A. Richter, M. Pleimling and A. Hiiller, cond-mat/0503733 
[17] H. Behringer, J. Phys. A: Math. Gen. 37, UA'i (!2004). 

[18] H. Behringer, On the structure of the entropy surface of microcanonical systems (Mensch und 

Buch Verlag, Berlin, 2004). 
[19] M. D'Agostino, F. Gulminelli, Ph. Chomaz, M. Bruno, F. Cannata, R. Bougault, F. Gramegna, 

I. Iori, N. Le Neindre, G. V. Margagliotti, A. Moroni and G. Vannini, Phys. Lett. B 473, 219 

(2000). 

[20] M. Schmidt, R. Kusche, T. Hippler, J. Donges, W. Kronmuller, B. von IssendorfT, and H. 

Haberland, Phys. Rev. Lett. 86, 1191 (2001). 
[21] H. Behringer, J. Phys. A: Math. Gen. 36, 8739 (2003). 
[22] K. Binder, Physica 62, 508 (1972). 
[23] D. P. Landau, Phys. Rev. B 13, 2997 (1976). 

[24] M. N. Barber, in Phase Transitions and Critical Phenomena, vol 8, edited by C. Domb and 
J. L. Lebowitz (Academic Press, London and New York, 1983). 



